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ABSTRACT 

We present results of a study aimed at deriving fundamental properties of cir- 
cumstellar debris disks from observed infrared to submillimeter spectral energy 
distributions. This investigation is motivated by increasing telescope/detector 
sensitivity, in particular the expected availability of the Space Infrared Telescope 
Facility (SIRTF) followed by the Stratospheric Observatory for Infrared Astron- 
omy (SOFIA), which will enable detailed studies with large source samples of late 
stage circumstellar disk and planetary system evolution. We base our study on an 
analytic model of the disk density distribution and geometry, taking into account 
existing constraints from observations and results of theoretical investigations of 
debris disks. We also outline the effects of the most profound characteristics of 
circumstellar dust including the grain size distribution and dust chemical com- 
position. In particular we find that an increasing iron content in silicates mainly 
causes an increase of the dust absorption effiency and thus increases the dust 
reemission continuum. Furthermore, the influence of the sp 2 /sp 3 hybridization 
ratio in carbon grains on the spectral energy distribution is examined. We in- 
vestigate the influence of various parameters on the resulting dust scattering and 
absorption/reemission spectral energy distribution and discuss the possibility for 
distinguishing between different disks from their infrared to submillimeter spec- 
tra. The strength and shape of amorphous silicate may be particularly diagnostic 
of debris disk evolutionary stages. Since the appearance of these features at 10/im 
and 20/im depends on the relative abundance of small grains - and therefore the 
minimum grain size and slope of the grain size distribution - they can be used 
to trace recent collisional processes in debris disks. Thus, debris disk surveys 
containing statistically large numbers of objects should reveal the likelihood of 
collisions and therefore the evolution of dust/planetesimals in debris disks. The 
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results of our study underline the importance of knowledge of the stellar photo- 
spheric flux, especially in the near to mid-infrared wavelength range, for a proper 
analysis of debris disk spectral energy distributions: While the quality of subtrac- 
tion of the direct stellar light at far-infrared wavelengths determines the accuracy 
of the mass estimate in the disk, our simulations show that the remaining stellar 
contribution due to scattering at near- to mid-infrared wavelengths constrains 
the dust grain size and chemical composition, e.g. the iron abundance in silicate 
grains. 

Subject headings: radiative transfer — methods: numerical — (stars:) circum- 
stellar matter — (stars:) planetary systems: protoplanetary disks 

1. Introduction 

Debris disks are solar system sized dust disks with micron-sized grains produced as 
by-products of collisions between asteroid-like bodies left over from the planet formation 
process. In the case of our solar system, the debris of Jupiter-family short-period comets 
and colliding asteroids represents the dominant source of zodiacal dust located between Mars 
and Jupiter. A second belt of dust is located beyond the orbit of Neptune (see, e.g., Dermott 
et al. 1992, Liou, Dermott, & Xu 1995). Besides the solar system, optical to mid-infrared 
images of (3 Pic (see, e.g., Kalas & Jewitt 1995; Weinberger, Becklin, & Zuckerman 2003) and 
submillimeter images of Vega, Fomalhaut, and e Eri (Holland et al. 1998; Greaves et al. 1998) 
have revealed spatially resolved debris disks which were first inferred from observations of 
infrared flux excesses above photospheric values with IRAS. Based on studies with ISO, the 
disk fraction is thought to decrease dramatically with age, amounting to much less than 
10% for stars with ages > lGyr (e.g., Spangler et al. 2001; see also Habing et al. 2001). 
Planetary debris disks are assumed to represent the almost final stage of the circumstellar 
disk evolution process, i.e., they are the evolutionary products of ongoing planet formation. 

In contrast to optically thick young circumstellar disks around Herbig Ae/Be and T Tauri 
stars with spatial structures dominated by gas dynamics, the much lower optical depth and 
lower gas-to-dust mass ratio in debris disks (Zuckermann, Forveille, & Kastner 1995; Dent et 
al. 1995; Artymowicz 1997; Liseau & Artymowicz 1998; Greaves, Coulson, & Holland 2000; 
Thi et al. 2001) let the stellar radiation - in addition to gravitation - be responsible for the 
disk structure. The Poynting- Robertson effect, stellar wind drag, and - if existing - gravi- 
tational stirring by embedded planets are all important in determining the dust population 
and disk structure (Liou & Zook 1999, Grady et al. 2000, Moro-Martm & Malhotra 2002). 
Similar to T Tauri-like disks, however, embedded planets may alter the debris disk structure 
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remarkably (see, e.g., Gor'kavyi et al. 1997, Kenyon & Bromley 2001). While planets may 
open gaps in gas-dominated young circumstellar disks (see, e.g., Bryden et al. 1999, 2000; 
Kley 2000), they decrease the particle density within their orbits in the case of a debris disk, 
provided that the dust sources are located outside the planet's orbit. 

Since the mass of small grains in debris disks and therefore the thermal dust reemission 
from these disks is much smaller than in case of T Tauri disks, only a very limited sample 
of observations exists so far. However, because of the high sensitivity of the mid-infrared 
detectors aboard the Space Infrared Telescope Facility (SIRTF) a substantial increase in the 
total number and in the specific information about debris disks is expected (c.f. Meyer 2001). 
In order to allow a quick classification of the basic characteristics of observed objects we 
perform a preparatory study with the goal of revealing the influence of major disk and 
dust parameters, such as the inner/outer radius, density distribution, grain size distribution, 
composition, etc., on the emergent spectral energy distributions (SED) of debris disk systems. 

Our first investigation, described in this publication, is based on a simply-structured, 
optically thin analytic disk model with a limited number of basic parameters. This approach 
allows us to distinguish clearly between the influences of different model parameters. Fur- 
thermore, it provides an easy to reproduce, flux-scalable (flux is directly proportional to dust 
mass in the very optically thin limit) database of debris disk SEDs 1 . In order to consider a 
reasonable model setup and parameter ranges we use boundary constraints given by previous 
observations and numerical simulations. 

In § 2 we describe the basic model and the method by which the SEDs have been cal- 
culated. In § 3, § 4, and § 5 we discuss specific properties of different species expected in 
debris disks and provide a detailed analysis of the influence of the different model parame- 
ters on the SED. Finally, we derive the SIRTF-appropriate colors m(8/im) — m(24/im) and 
m(24/im) — m(70//m) from the simulated debris disk SEDs (§ 6). 

2. Model description 

According to the expected structure of debris disks, we base our simulations on the 
assumption of an optically thin disk, whereby the optical depth is measured in the disk 
midplane as seen from the star. In this case, multiple scattering of radiation and dust heating 
due to dust-reemission can be neglected. The only radiative processes to be considered 



1 A database with all dust SEDs presented in this publication is available at 
http://mc.caltech.edu/~swolf/downloads/debris/ (or contact the authors) 
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are scattering, absorption, and reemission of stellar radiation by dust grains. The two- 
dimensional radiative transfer problem can then be reduced to a one-dimensional problem. 

In all simulations we chose the shortest wavelength of stellar emission to measure the 
optical depth. This wavelength amounts to 0.2 jum in practice due to the focus on solar-type 
stars and the availability of optical constants. In order to satisfy the optically thin assump- 
tion in our calculations we apply a limit for the optical depth of r max = 0.05, which translates 
into a maximum attenuation of stellar light by 4.8 % at 0.2 /im. Depending on the grain size 
and chemical composition, the absorption efficiency decreases towards longer wavelengths, 
but may reach a local maximum in the mid-infrared wavelength region for particular dust 
chemistries (see, e.g., Fig. 3 and Fig. 5 for the wavelength-dependent absorption efficiency 
of MgSi0 3 and crystalline Olivine, respectively). However, the flux of the considered stel- 
lar radiation source is small in the mid- infrared compared to that in the ultraviolet /optical 
wavelength region (S„(10fim)/S^(0.55fj,m) 1.3 x 10~ 4 for the solar SED: see below). Due 
to the small amount of dust in the disk/shell for r max = 0.05, the flux reemitted in the 
mid-infrared wavelength range by hot dust is also negligible for further disk heating. 

Each dust grain is heated by direct stellar radiation only. Thus the dust grain tem- 
perature is a function of the optical parameters of the grains, the incident stellar radiation, 
and the distance d from the star. In this case the radiative transfer equation has a simple 
solution which allows one to derive the distance from the star at which the dust has a certain 
temperature. Let 

L A , S = 47ri2 B 2 7rS A (7;) (1) 

be the monochromatic luminosity of the star (radius R s , effective temperature T s ) at wave- 
length A and 

= 4na 2 Q? s nB x (T g ) (3) 

be the absorbed and reemitted luminosity of a dust grain with radius a and resulting tem- 
perature T g at the (unknown) distance d from the star. Using the constraint of energy 
conservation 
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one derives the distance of the grain from the star as 

(5) 



d(T g ) = f 
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If the dust sublimation temperature is known, Eq. 4 also allows one to estimate the subli- 
mation radius for each dust component in the shell (characterized by the grain radius and 
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chemical composition). The flux of the light scattered by a single dust grain amounts to 



where A is the dust grain's albedo. 

We calculate the temperature distribution of the dust based on 500 logarithmically 
equidistantly distributed wavelengths in the interval [0.2/xm, 500/xm]. According to Eq. 5, 
distances d(T) corresponding to 500 temperatures equidistantly distributed between 2.73 K 
and the dust grain sublimation temperature of the grains are determined. The dust reemis- 
sion and scattering is calculated at 200 logarithmically equidistantly distributed wavelengths 
in the interval [2/iin, 200/xm]. The latter interval was chosen according to the wavelength 
range covered by broad band filters aboard SIRTF: 3.6/im - 160//m. The net spectral energy 
distribution results from a simple summation of the reemitted and scattered light contribu- 
tions from all grains. 

Since the SED depends on the radial density distribution only, we consider a spheri- 
cal shell instead of a circumstellar disk. This simplification has no effect on the resulting 
temperature distribution and/or SED, but reduces the number of free parameters. If a disk 
with an opening angle from the midplane 9 - but the same radial density distribution - is 
considered instead, the SEDs shown in the subsequent sections have to be multiplied by the 
factor 



but the relative shape of the dust contribution to the net SED by scattering, absorption and 
reemission remains unchanged. For 9 = 10°, f$ m 0.087 while for 9 = 45°, f$ m 0.38. 

Our parameters are as follows, unless specified otherwise in subsequent sections. The 
applied radial density profile in the disk/shell follows a power-law 



Except for a parameter study in § 4.1 we use q — 1, representing a debris-disk with no 
perturbations by embedded planets and no particular assumptions about dust production 
processes (see § 4.1 for a more detailed description). If not otherwise noted, the mass of 
the disk amounts to 10~ 10 M Q , the inner disk radius is equivalent to the dust sublimation 
radius, and the outer radius amounts to 100 AU. The central star has a solar-type SED (we 
use the solar SED published by Labs & Neckel 1968, extended by a blackbody SED beyond 
151 jitm). The distance to the system is assumed to be 50 pc. 




(6) 
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3. Dust grain parameters 

The dust found in young, massive circumstellar disks around TTauri and Herbig Ae/Be 
stars consists mainly of two distinct species: a silicate component and a carbonaceous com- 
ponent (Savage & Mathis 1979, Draine & Lee 1984; Malfait et al. 1999). Based on the 
investigation of the interstellar extinction Weingartner & Draine (2001; see also Draine & 
Lee 1984) derived optical properties and relative abundances of both. Despite a different 
size distribution of the grains in the circumstellar vs. the interstellar environment (see, e.g., 
Mathis, Rumpl, & Nordsieck 1977; Beckwith & Sargent 1991; Miyake & Nakagawa 1993), 
the chemical composition turned out to be very similar, allowing detailed studies of the 
circumstellar environment and disks of young stellar objects. However, the very complex, 
(magneto-) hydrodynamically dominated optically thick disks permits deriviation of only 
very basic further information about the dust, such as the grain size distribution exponent 
and minimum and approximate maximum upper grain size (e.g., Wolf, Padgett, & Stapelfeldt 
and references therein). In contrast to the situation for young optically thick disks, older 
optically thin debris disks are expected to allow a more detailed investigation of dust grain 
parameters, in particular the chemical components and the grain sizes. An important caveat 
to this claim of simpler modelling for debris disks is that because the dust dynamics are 
determined by a combination of radiation pressure, the Poynting-Robertson effect, and the 
(possible) production of dust due to collision processes and gravitational interaction with 
possibly embedded planets, the disk density distribution may be dictated by several addi- 
tional parameters, the influence of which on the dust reemission/scattered light SED will be 
discussed in detail in Sect. 4 and explored more thoroughly in a forthcoming paper. 

Fortunately, existing laboratory measurements provide a rich database of dust properties 
expected in different environments (defined by the dust formation conditions, temperature, 
abundances of particular elements, etc.). A compilation of a large number of measured optical 
dust parameters has been published by Henning et al. (1999). Before applying selected 
chemical compositions in our investigation of the influence of disk model parameters on the 
observable SED, we outline in brief the most characteristic differences between the dust 
compositions expected to dominate. 

We consider the following chemical components: 

1. Amorphous silicates with varying iron content from Fe-poor through Fe-rich, 

2. Crystalline silicates (in particular Olivine), and 

3. Amorphous carbon species with different graphite/diamond structure ratios. 
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ISO spectroscopy of Herbig Ae/Be stars revealed several further abundant chemical com- 
ponents, such as Fe (broad, weak emission feature at ~ 2 — 4/xm; see Henning & Stog- 
nienko 1996), FeO (broad emission features at ~ 21 — 25/im; see Henning et al. 1995) 
and/or FeS (see Henning & Stognienko 1996) and H 2 ice (broad emission features be- 
tween pa 40 — 80/im; Warren 1984). For model SEDs taking into account these additional 
components see, e.g., Malfait et al. (1999), Bouwman et al. (2000), and Meeus et al. (2001). 
However, the SEDs of debris disks measured so far do not allow one to perform a comparably 
detailed chemical analysis; SIRTF is expected to be revolutionary in this regard. 

We consider the dust grains to be homogeneous spheres. Although dust grains are 
expected to have a fractal structure, the scattering behaviour is similar to that of spheres 
(see Lumme & Rahola 1994 for porous dust particle light scattering). Furthermore, dust 
grains are expected to have a non-spherical shape (see, e.g., Elvius & Hall 1967; Scarrott, 
Draper, & Warren-Smith 1989; Hajjar & Bastien 1996; Kastner & Weintraub 1996; Dollfus 
& Suchail 1987; Johnson & Jones 1991; Chrysostomou et al. 2000). However, we do not 
consider the grains to be aligned on a large scale by magnetic fields and thus, the assumption 
of spherical grains is a valid approximation (Wolf, Voshchinnikov, & Henning 2002). The 
optical parameters of the grains required for the estimation of the dust absorption, reemission 
and scattering, are derived from the complex refractive index m = n + ik and the grain 
radius. We use measurements of the complex refractive index m published by Dorschner 
et al. (1995; Silicates) and Jager et al. (1997; Carbon) 2 . Furthermore, we compare selected 
dust properties with the optical properties of "astronomical silicate" and graphite provided 
by Weingartner & Draine (2001; see also Draine & Lee 1984) since these data are used in 
a large variety of simulations of circumstellar dust configurations. According to Eqs. 5 and 
6, the absorption and scattering efficiencies determine the particular contribution of each 
chemical component to the net SED. We consider absorption/re-emission and scattering for 
individual grains of given composition and size over the wavelength range 0.2 - 200/zm. 

The interaction of the stellar radiation field with the dust grains is described by Mie 
scattering theory. We calculate the Mie scattering function using the numerical solution for 
the estimation of the Mie scattering coefficients published by Voshchinnikov (2003), which 
achieves accurate results both in the small as well as in the - arbitrarily - large size parameter 
regime 3 



2 The complex refractive indices are available at 
http : //www . astro . uni- j ena . de/Laboratory/Database/odata . html. 

3 The code for calculation of the Mie scattering coefficients is available at 
http : //mc . caltech . edu/^swolf /miex-web/miex . htm . 
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In the following we briefly describe the silicate and carbon species considered: 

Silicates: Previous studies concerning the mid-infrared SED of evolved stars and the 
(3 Pic circumstellar disk revealed that a large variety of silicate species is present in these 
different environments. For example, Pantin, Lagage, & Artymowicz (1997) find amorphous 
olivine (MgFeSi0 4 ) to be the dominant chemical dust component in the (3 Pic disk (55% 
amorphous olivine, 35% amorphous pyroxene, and 10% crystalline Olivine). In contrast 
to this, Molster, Waters, & Tielens (2002a,b) find narrow-band emission features in the 
SED of dust around evolved stars, clearly pointing to the existence of crystalline silicates. 
In order to cover this large range of different silicates, we investigate the main differences 
between a number of different amorphous pyroxene species (Mg^Fei-^SiOs with £=0.4, 0.5, 
0.6, 0.7, 0.8, 0.95 and 1.0), and amorphous Olivine species (Mg2j / Fe 2 -2 2 /Si04 with £=0.4 and 
0.5). The efficiency factors for Ifim. grains are shown in Fig. I. Despite the characteristic 
10/iin and 19/im maxima of the absorption efficiency, the most significant property is the 
strong dependence of the ultraviolet (UV) to near-infrared (NIR) absorption efficiency on 
the Fe content (see Jager et al. 1994 for a detailed compilation and discussion of absorption 
bands of amorphous pyroxene). With increasing relative Fe content the absorption efficiency 
increases by about three orders of magnitude in the case of a micron-sized grain in this 
wavelength range from Fe-deficient MgSi0 3 to Fe-rich Mg 8 Fei. 2 SiO4. Thus, the temperature 
and reemitted flux of the grains at a given radial distance from the star increases. These 
effects are demonstrated in Fig. 2 where the resulting density distribution for our standard 
disk model (see below) is shown. Furthermore, the comparison with "astronomical silicate" 
(Weingartner & Draine 2001) shows that it is much more similar to Fe-rich than to Fe-poor 
silicates. For comparison, the scattering/reemission SED of crystalline Olivine (MgFeSiO^ 
is shown. Similar to Fe-poor silicates, crystalline Olivine has a very low absorption efncieny 
in the wavelength range of stellar emission (see Fig. 5) and can therefore only be traced 
(by its pronounced numerous narrow emission features) either based on highly accurate 
photometric measurements or if a high relative amount of this material is present in the 
disk. Furthermore, the dependence of the efficieny factors on the grain size is illustrated in 
Fig. 3, 4, and 5 for Fe-deficient and Fe-rich amorphous silicate and crystalline Olivine. 

Carbon: It is commonly assumed that the carbonaceous dust component occurs in 
the form of graphite (see, e.g., Draine & Lee 1984, Li & Greenberg 1997, Weingartner & 
Draine 2001). In contrast to this common assumption, the dust emissivity described by 
a power-law with a spectral index of « 1 (Campbell et al. 1976; Sopka et al. 1985; Mar- 
tin & Rogers 1987; Giirtler, Kompe, & Henning 1996) can be explained by the occurence 
of very disordered two-dimensional material like amorphous carbon (Kittel 1963). Jager, 
Mutschke, & Henning (1998) synthesized structurally different carbon material with an in- 
creasing sp 2 /sp 3 ratio by pyrolizing cellulose materials at 400-1000°C The different carbon 
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structures therefore represent different graphite (sp 2 hybridization) to diamond structure 
(sp 3 hybridization) abundance ratios. In Fig. 6 and 7 we show the efficiency factors calcu- 
lated for the two extreme cases: the 400 K carbon modification (graphite-poor) and 1000 K 
carbon modification (graphite-rich). As these figures show, carbon has a high absorption 
efficiency in the UV/visual wavelength range, but - in contrast to silicate species - no promi- 
nent absorption features (see also Fig. 2). Thus, it mainly adds a strong continuum to the 
SED. The absorption efficiency in the UV/visual wavelength range is comparable to Fe-rich 
silicates. With increasing sublimation temperature the range of almost constant, high ab- 
sorption efficiency is extended towards longer wavelengths. The particular turn-over point 
beyond which the absorption drops continiously, increases as the grain size is increased. 
It is noteworthy that the absorption efficiency of graphite as published by Weingartner & 
Draine (2001) is smaller than that of the other carbon configurations in the mid-infrared 
wavelength range. 

Finally, we illustrate the importance of scattering of stellar radiation by the dust in the 
case of iron-deficient and iron-rich silicate (see Fig. 8). Mainly determined by the solar-type 
SED of the star, the contribution of scattering to the resulting SED is negligible beyond 
w 10/xm. At shorter wavelengths, however, the lower absorption efficiency in the case of an 
iron-deficient silicate results in a larger relative constribution of the scattered light, while 
the iron-rich amorphous Olivine dust reemission dominates at wavelengths of a few microns. 
Thus, if iron-poor silicates or crystalline Olivine are the dominating dust species, spectral 
features of the star will be present in the near- to mid-infrared wavelength range even after 
correct substraction of the stellar SED (e.g., resulting from independent stellar models). 

4. Parameter Study 

In this section we explore the parameter space allowed by our adopted analytic disk 
model. We are strongly guided by observations and theoretical constraints in considering 
viable values for the disk density distribution, inner and outer disk radius, and the grain size 
distribution. 



4.1. Disk density distribution 

Based on hydrostatic models, the radial density distribution in T Tauri disks is proposed 
to be in the range 

n(r) oc r -( L9 - 2 - 4 ) (9) 
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(see, e.g., Kenyon & Hartmann 1987, D'Alessio et al. 1999, Chiang & Goldreich 1997). De- 
tailed modelling of the near-infrared - to - millimeter appearance of several spatially resolved 
TTauri disks has confirmed these predictions (see, e.g., Burrows et al. 1996; Stapelfeldt et 
al. 1998; Cotera et al. 2001; Wood et al. 2002; Wolf et al. 2003). Optically thin debris disks 
which - in contrast to TTauri disks - are not dominated by gas, are assumed to cover this 
range of radial density profiles as well but also to manifest values outside of this range. 

In the most simple model of debris disks, structure is dominated exclusively by grav- 
itation, radiation pressure, and Poynting-Robertson drag. In a collisionless system with- 
out sources or sinks and grains in circular orbits, the exponent q amounts to 1 (see, e.g., 
Briggs 1962). This "classical solution" approximately represents the overall distribution of 
dust in the solar system (see, e.g., Ishimoto 2000). However, as for instance Gor'kavyi et 
al. (1997) pointed out, there exist several effects that may change this distribution remark- 
ably, including resonance effects with planets, gravitational encounters with planets which 
occur in the form of elastic gravitational scattering, mutual collisions of particles, evapora- 
tion of dust grains, and existence of sources of dust with highly eccentric orbits (such as the 
Encke comet in the case of the Solar system; (see, e.g., Whipple 1976, Sykes 1988, Epifani 
et al. 2001). The resulting density distribution is in most cases no longer described by a 
single power-law as stated in Eq. 8, but depends strongly on the properties (mainly the or- 
bits and masses of planets and dust sources) of each particular disk/planetary system. The 
exponents of the radial density distribution derived for different radial zones in the model of 
the Solar system reach from g=1.0 to 2.4 (Gor'kavyi et al. 1997). Another simply-structured 
debris disk would be a cloud of grains moving outward from the centre in hyperbolic trajec- 
tories. This solution, which applies to so-called /3-meteorites in the interplanetary space is 
described by q=2 (see, e.g., Lecalvier des Etangs, Vidal-Madjar, & Ferlet 1998; Ishimoto & 
Mann 1999). In the case of (3 Pic, visible observations of the scattered starlight and mid- 
infrared images show that the outer disk (at radii > 100AU) can be described by a power-law 
in the range [1.5,2.3] (Artymowicz, Burrows, & Paresce 1989; Kalas & Jewitt 1995; Smith 
& Terrile 1984; Pantin et al. 1997) 4 . In the following we consider the cases q—1.0, 1.5, 2.0, 
and 2.5. 

We find that for a given grain size and chemistry, the radial density profile, described 



4 We want to remark that Trilling et al. (2000; see also Trilling & Brown 1998) reported the detection of 
circumstellar disks around 55 Cnr, p CrB, and HD 210277 which are planet-harbouring systems (Butler et 
al. 1997, Noyes et al. 1997, Marcy et al. 1999). From scattered-light images these authors derive a volume 
density profile described by q = 4. However, since these measurements could not be confirmed so far by later, 
independent (see, e.g., Jayawardhana et al. 2002) measurements, we do not take so steep density profiles 
into account. 



- 11 - 



by the exponent q (see Eq. 8), determines the overall slope of the SED (see Fig. 9). The 
relative strength of the spectral features is not affected by q. The steeper the density profile, 
i.e., the larger the exponent q, the larger the relative amount of dust located in the inner 
regions of the disk. Since the dust temperature increases towards the star, this results in an 
increase of the near/mid-infrared-to-(sub)mm luminosity ratio. But also the net reemission 
of the disk increases with increasing steepness of the density profile since the relative number 
of hotter and therefore more luminous grains increases. 

4.2. Inner Gaps in the Disk 

Observed SEDs of selected disks around T Tauri and debris-type disks have been found 
to show hints of inner cavities such that the inner radius of the disk is apparently much 
larger than the sublimation radius of possible dust species. Examples among T Tauri disks 
are GM Aurigae (Koerner et al. 1993, Rice et al. 2003) and TW Hya (Calvet et al. 2002). In 
several prominent debris disks, inner cavities have been found as well: f3 Pic (inner radius: 
20 AU 5 ), HR 4796A (30-50 AU), e Eri (50 AU), Vega (80 AU), and Fomalhaut (125 AU) - see, 
e.g., Dent et al. 2000; Greaves, Mannings, & Holland 2000b; Wilner et al. 2002; Holland et 
al. 2003. Taking into account the physical processes responsible for the spatial dust density 
distribution in debris disks, scattering of particles by massive planets is assumed to be the 
major effect in explaining these large inner holes: Dust grains drifting inwards due to the 
Poynting- Robertson effect are likely to be scattered into larger orbits resulting in a lower dust 
number density within the planet's orbit. Based on the assumed cavity sizes of the debris 
disks quoted above (relative to the total sizes of these disks), we consider the following cases: 
0.1 AU, 1.0 AU, and 10 AU and compare the resulting SED to the case of a disk without any 
gap, that is, extending inward to the dust sublimation radius. 

Clearing the inner disk results in a loss of warm dust which is mainly responsible for 
the near/mid-infrared shape of the SED. With increasing gap size the minimum value of the 
mid- infrared flux becomes smaller and is shifted towards longer wavelengths (see Fig. 10). Be- 
cause of the higher temperature and therefore more significant contribution to the near/mid- 
infrared spectrum by small grains, the decrease of the flux in this wavelength range is more 
pronounced when the relative number of small grains in the system goes up. Keeping the 
disk mass constant, the flux in the millimeter region (not accessible by SIRTF) increases 



5 We want to remark that this value derived by Dent et al. (2000) from optical-to-submillimeter spec- 
troscopy is in contrast to the result given by Weinberger et al. (2003) who assume a substantial amount of 
dust grains within that radius based on spatially resolved mid-infrared spectroscopy (see § 4.4). 
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slightly, but the net flux is smaller compared to a disk without an inner gap. This is because 
the fraction of the stellar flux absorbed by a single dust grain decreases with increasing radial 
distance from the star (see Eq. 2). Applying the same explanation to the scattering cross 
section (see Eq. 6), the slight decrease of the scattered stellar light with increasing gap size 
can be explained by analogy. 

4.3. Outer disk radius 

Similar to T Tauri disks, circumstellar disks at the late stage of their evolution appear to 
show a large range of possible radii, typically in the range of about 100 AU to several hundered 
AU. So far, scattered light images and SEDs revealed radii of 80 AU (eEri), 120 AU (Vega), 
125 AU (HR 4796), 185 AU (Fomalhaut) and 1000 AU in case of the disk around (3 Pic (see 
Vidal-Madjar et al. 1998; Dent et al. 2000; Greaves, Mannings, & Holland 2000b; Holland et 
al. 2003). However, scattered light images as well as near /mid- infrared SEDs of debris disks 
are mainly determined by the pronounced spectral features and the total flux in the small 
grain component in the disk (see § 4.1, Fig. 12), while larger grains with radii of several 
tens of micron and larger can be more efficiently found by (sub)millimeter measurements. 
One has to take into account that the abundance and spatial distribution of the small grain 
component, especially the one with f3 > 1 (see Eq. 11), depends on the location of dust 
sources and the likelihood of collisional events. Since these parameters are not necessarily 
correlated with the actual disk size, sub millimeter/millimeter measurements are required to 
trace larger grains which are less affected by the radiation pressure and Poynting-Robertson 
effect. The question to be answered here is whether there are significant differences in the 
SED resulting from different disk sizes. 

As result of our simulations, Fig. 1 1 shows the SED for disks with outer radii between 
25 and 400 AU. The influence on the near /mid- infrared wavelength range is seen only as a 
decrease of the net flux with increasing disk size which is due to the decrease of the mean 
temperature of the disk. At several tens of microns and beyond lOOyum also the shape of 
the SED is affected because the smaller mean temperature results in an increase of the 
(sub) millimeter-to-near /mid-infrared flux ratio. 

4.4. Grain size distribution 

In this section we consider dust models with a grain size distribution which is charac- 
terized by a (minimal) set of three parameters: the minimum and maximum grain radius 
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(a m i n , a max ) and a power-law exponent p of the size distribution: 

n(a) oc a~ p . (10) 

Since debris disks are assumed to be a place of continuous dust production due to collisional 
processes between larger particles (rather than being nurtured from a hostile interstellar 
medium - see, e.g., Artymowicz & Clampin 1997), all three parameters depend on the par- 
ticular debris disk system. In contrast to the interstellar medium or disks around very young 
stellar objects, all three parameters may therefore vary in a wide range even between different 
radial zones within one and the same disk (Gor'kavyi et al. 1997). Furthermore, the collision 
rate is assumed to decrease as the disk ages since the fraction of disk material present in 
larger bodies increases with time (see, e.g., Kokubo & Ida 2002). This might allow us to 
classify debris disks based on their grain size distribution, at least on a statistical basis. 

Constraints on the minimum grain size a min : The smallest grain size expected to be present 
is determined on the one hand by the likelihood of collisional processes in the disk, and on 
the other hand by the optical properties of the dust and the strength of the stellar radiation 
field. If no collisions are occuring, the minimum grain size can be estimated from the ratio 
of the radiation pressure to gravity, also known as (3, (Burns et al. 1979) 

0=— LsQpT (11) 

16?r GM s cap g ' v ; 

where L s /M s is the central stellar luminosity/mass, G is the gravitational constant, c is the 
speed of light, p g is the material density of the grain, and Q pr is the radiation pressure 
coefficient averaged over the entire stellar spectrum. Those grains for which (3 > 1 are 
"blown away" by radiation pressure on a timescale much shorter than the age of the star 
(e.g., Aumann et al. 1984). Mid-infrared observations of the (3 Pic disk by Weinberger et 
al. (2003) show that there must exist a substantial amount of grains with radii < 10pm within 
20 AU, while Pantin et al. (1997) even constrain the minimum grain size to be as small as 
0.1/im. In case of Fomalhaut's debris disk, Wyatt & Dent (2002) derive a minimum grain 
size of ~ 7pm from submillimeter observations. In our parameter study we consider a min 
= 0.1pm and 10/im since having the bulk of the mass in larger grains does not produce a 
substantial mid-infrared contribution to the SED. 

Constraints on the maximum grain size a max : The upper grain size is a very uncertain param- 
eter since debris disks are expected to harbour dust "particles" ranging from (sub) micron- 
sized grains to planetesimals/planets as the result of an (almost) finished process of dust 
grain growth, planetesimal formation, and simultaneous dust production due to collision 
events (see, e.g., Beckwith, Henning, & Nakagawa 2000). However, according to the wave- 
length range considered, grains larger than « 1mm will not give a remarkable contribution 
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to the SED. The existence of these grains can therefore not be proven or ruled out by SIRTF 
measurements. Complementary ground-based (sub) millimeter observations can in principal 
increase this upper grain size limit by about one order of magnitude. Based on previous 
observations of the very limited number of debris disks only very weak constraints on upper 
grain sizes could be derived. Paresce & Burrows (1987), for example, found an absence of sig- 
nificant colour excess in visible observations at the outer part of the /3 Pic disk, which implies 
that the dominant particles for the scattering are larger than a few /xm. Furthermore, ob- 
servations of main-sequence and pre-main-sequence (post-Herbig Ae/Be, post-T Tauri) stars 
show that the distribution is broad and extends from below 1 /im to at least ~1 mm (see, 
e.g., Skinner, Barlow, & Justtanont 1992; Sylvester et al. 1996; Sylvester & Skinner 1996; Li 
& Greenberg 1998). In our parameter study we consider a max = 0.1mm, 1mm, and 10mm. 

Constraints on the grain size distribution exponent p: The equilibrium size distribution re- 
sulting from a collisional cascade can be described by a power-law (Eq. 10) with the exponent 
P = 3pd _ 2 (Dohnanyi 1969). Furthermore, Tanaka, Inaba, & Nakazawa 1996 could show 
that pd = 1.833 (p = 3.499) describes the case of a self-similar (and therefore infinite) 
collisional cascade. Wyatt & Dent (2002) find in fact a very good agreement of the dust 
SED from Fomalhaut's debris disk with the assumption of p = 3.499 (see also Holland et 
al. 2003). However, in general, real dust particle/planetesimal ensembles are expected to 
be characterized by size-dependent particle strengths, therefore resulting in a distribution of 
exponents p as a function of grain size. The "mean" size distribution exponent could then 
be slightly shifted from 3.499 (Durda & Dermott 1997). For example, Patin et al. (1997) 
derived two different exponents to fit mid-infrared images of the f3 Pic disk: p=3.0 and 
3.3 for the smallest grains (with radii up to 10//m) and the larger population, respectively. 
Since particle growth/collision processes are expected to be mainly dependent on the radial 
distance from the star (Weidenschilling & Davis 2000), the dependence on the (mean) grain 
size translates into a dependence of the exponent p on the location within the disk. Another 
argument for a clear deviation from p = 3.499 comes from the fact that newly created small 
grains are likely to be removed from the system faster than larger ones may decrease this 
value, at least on short timescales, i.e., as long as a small grain population resulting from 
collision processes is present (see, e.g., Moro- Martin & Malhotra 2003). We consider two 
different values of the exponent, p = 2.5 and 3.5. 

In Fig. 12 the SEDs resulting from dust disks with different grain size distributions are 
shown. Based on the discussion above, we consider the following parameter combinations: 
a min = 0.1/xm, 10/im; a max = 0.1mm, 1mm, 10mm; and p = 2.5 and 3.5. Comparing the 
SEDs within these parameter ranges and the considered wavelength range one can draw the 
main conclusion that the relative number of small grains in the system (determined by a m i n , 
a max , and p) influences the shape of the SED. In particular we find that: 
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• Increasing the maximum grain size results in a decrease of the net flux and a decrease 
of dust composition specific emission features. Since the amount of the underlying 
continuum can be due either to almost featureless emission by carbon grains of arbitrary 
size (distribution) or large silicate grains, the maxium grain size might become the 
most difficult parameter to be derived from SIRTF measurements alone. However, 
millimeter measurements can be used to distinguish between the scenarios based on 
different slopes of the SED (similar to what has been done in case of T Tauri disks; see 
Beckwith & Sargent 1991). 

• Given a fixed disk mass, increasing the grain size distribution exponent and decreasing 
the minimum grain size increases the relative amount of small grains and therefore 
pronounced dust-specific spectral features in the SED. 

The SEDs shown in Fig. 12 illustrate that the relative abundance of small grains and 
therefore the minimum grain size and grain size distribution exponent will be well-constrained 
disk parameters from SIRTF measurements due to the abundance of spectral features in the 
5-40/im range (see also Fig. 3, 4, 5). 

5. SED probability distribution 

In this section we assess the likelihood of observing a particular SED given an assumed 
underlying distribution of the parameters we have varied in the previous subsections. We 
consider a plausible range of values for the minimum and maximum grain size, the radial 
density distribution exponent q, and the grain size distribution exponent p. The other disk 
and dust parameters are fixed: the inner radius (= sublimation radius), the outer radius (= 
100 AU) and the chemical composition of the dust (amorphous Olivine - MgFeSi0 4 ). Ac- 
cording to the previous parameter range justifications for the variable quantities we calculate 
first the separate SEDs for 10 3 grain sizes logarithmically distributed between 0.1/zm and 
1000/iin in disks with 31 different radial density distribution exponents linearly distributed 
in the range [1.0,2.5], i.e., 3.1 x 10 4 separate SEDs. In the second step we select a lower 
and upper grain size randomly from this grain size distribution, assuming an equipartition 
of grain sizes on the prepared logarithmic grid of possible grain radii. Furthermore, we ran- 
domly chose the grain size distribution exponent p from one of 11 possible values linearly 
distributed in the range [2.5, 3.5] and one of the 31 radial density distribution exponents 
q. In the third step we derive the net SED for a disk defined by these randomly chosen 
parameter settings based on the linear weighted combination of the single SEDs calculated 
in the first step. In order to allow the comparison of the (shapes of) SEDs resulting from 
the 31 x 11 x X)i=i° % ~ 1.7 x 10 8 possible different disks, the net luminosity is normalized 
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to "1" (arbitrary unit). For this reason, the single SEDs had to be calculated over a larger 
parameter range than in the previous studies in order to minimize neglected flux reemission 
appearing outside the considered wavelength range; here we chose A = [0.2//m, 4000/xm]. In 
Fig. 13 the distribution of 1.2 x 10 7 normalized, individual disk SEDs is shown, whereby the 
density of SEDs in the flux-wavelength-space is color-coded. 

This SED density distribution of course depends on the particular distribution of dust 
and disk parameters within the given parameter ranges. In order to illustrate this, in Fig. 13 a 
second set of SED density distributions is shown for which the minimum and maximum grain 
size have been randomly selected from a linear distribution of grain sizes (instead of from a 
logarithmic one as before). The main difference between the approaches is the much stronger 
pronounciation of the silicate emission bands in case of a logarithmic distribution since in that 
case both the mean minimum and mean maximum selected grain size is smaller than in case of 
a linear grid of possible grain size radii. The "real" (unknown) distribution of possible grain 
size boundaries depends on the outcome of collision processes. The comparison of observed 
SED density distributions with distributions as shown in Fig. 13 but based on theoretically 
predicted collision scenarious may therefore provide a key to prove the underlying physics of 
collision processes in debris disks. 

The range in normalized flux density, however, over which possible debris disk SEDs are 
distributed does depend on the parameter range of the dust and disk parameters only, and 
is therefore the same in both apporaches. Our simulations show that in case of a possible 
grain size distribution in the range of 0.1-1000/zm the most narrow distribution of already 
normalized SEDs still spans a range over more than one order of magnitude in possible 
fluxes (in the range of A ~ 40 — 50/xm), but increases to above 3 orders of magnitude towards 
the lower and upper boundary of the considered wavelength interval. The scattering of the 
SEDs slightly decreases if the range of possible grain radii is decreased (also Fig. 13). In 
particular, we find an additional minimum of possible flux values around A ~ 50 — 70/im for 
our particular dust/disk setup in case of large grain size distributions (a = 10 — 1000/im). 
This spectral regime may be the most diagnostic for determining disk masses relatively 
independent of the details of the assumed disk or grain properties. 



6. SIRTF Colors 

Based on selected SEDs presented in previous sections we derive the colors m(8/im) — 
m(24/im) and m(24/xm) — m(70/xm), applying the filter characteristics of the SIRTF cameras 
IRAC (8//m) and MIPS (24/xm, 70/xm). The resulting color-color diagrams, shown in Fig. 15- 
14, are normalized to a blackbody with a temperature of 9650 K which sufficiently accurately 
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represents the mid-infrared - to - submillimeter SED of the AOV standard star Vega (Castelli 
& Kurucz 1994, Decin et al. 2003): 



m 1 -m 2 = -2.5 ( log 10 ^q^^'V - ) + (m h Standard - m 2 , Standard), (12) 



where 



mi, Standard ~ ?™2, Standard — — 2.5 log 10 ^O^ni „ ■ (13) 



g; m m ^^(9650K)rfA 

/oZr^geSO^A' 

Here, 5a = du/d\ is the dust scattering/reemission SED, Ti and T 2 describe the wavelength- 
dependent filter transmission 6 , and B\ is the Planck function. Note in these figures that only 
the dust emission / scattering spectrum is represented; the stellar photosphere has been re- 
moved in the SEDs. 

These diagrams reflect what has been outlined in the discussions about the possibility 
to distinguish between different disk models. For example, the color-color diagram for disks 
with different outer radii (Fig. 15, upper right plot, based on the SEDs given in Fig. 11) shows 
a large dynamic range of colors and therefore clearly distinguishable models only in case of 
a small dust grains (here: 0.1/zm). In contrast to this, a grain size distribution with grain 
radii of 10/im and 1mm the almost constant offset between the SEDs in the mid-infrared 
wavelength range results in a nearly constant value of m(8/im) — m(24yum). 



7. Summary and Conclusions 

In preparation for SIRTF and SOFIA observations of circumstellar debris disks we have 
investigated the influence of various disk and dust parameters on the resulting SED, in 
particular the offset from the stellar photospheric flux due to dust reemission and scattering 
processes. Based on an analytic disk model we considered parameter ranges as suggested from 
existing observations of a small number of debris disk systems as well as constraints derived 
from theoretical approaches. We restricted our investigations primarily to the wavelength 
range accessible by SIRTF. 

We draw the following main conclusions from our study of the parameter space: 

1. The correct estimation/substraction of the stellar photospheric flux is essential for 
the subsequent data analysis in two respects. First, it is important for the correct 



6 SIRTF instruments characteristics are available at http://sirtf.caltech.edu/SSC/obs/. 
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estimation of the dust continuum (not taking into account any particular dust emis- 
sion features) which mainly constrains the disk mass, the upper grain size and grain 
size distribution exponent (§ 4.4), the outer disk radius (§ 4.3), and the overall ra- 
dial disk density distribution exponent (§ 4.1). Secondly, the strength of remaining 
stellar spectral features allows one to derive conclusions about the scattering to ab- 
sorption efficiency (or albedo) of the dust grains and therefore - based on a simple 
silicate/carbon dust grain model - to derive constraints on the crystalline silicate / 
iron-poor amorphous silicate to iron-rich amorphous silicate / carbon ratio (§ 3). 

2. The characteristic amorphous silicate broad band solid state emission features are pre- 
dicted to be of decisive importance for the characterization of the evolutionary state 
of debris disks. The appearance of these features depends strongly on the existence 
of small grains (smaller than about 10/xm), the abundance of which, however, is as- 
sumed to be small in a collisionless system (depending on the stellar SED; § 3 and 
§ 4.4). The strength of these features will allow to decide if collision processes took 
place still recently in the disk of the particular star. However, one has to take into 
account that only "snapshots" of the current, in respect of the timescale of the colli- 
sional events arbitrarily selected states of debris disks can be obtained. The SED and 
thus the derived smallest size of highly abundant grains does therefore not necessarily 
describe the real, "mean" evolutionary state of a particular disk (if defined by the 
time-dependent collision rate). However, surveys based on a statistically large number 
of debris disk systems with large subsets of stars with a similar age will allow to derive 
general characteristics of debris disks as a function of their (evolutionary) age. 

3. The mid-infrared continuum will allow to conclude if an inner hole - due to the presence 
of at least one embedded planet - exists. Furthermore, the inner radius of that gap (= 
orbit of the planet) could be derived (§ 4.2). 



The presented study was performed to address the most fundamental questions one 
will want to answer with the analysis of the SED of debris disk systems. Based on this 
"first-level" analysis, a more detailed investigation of particular debris disk systems will be 
possible. For instance, the radial segregation of chemically different dust grains, the spa- 
tial distribution of ice-mantled grains with spectral features in the range of 10-100/im range 
(e.g., Warren 1984), the possibility to distinguish qualitatively different planetary architec- 
tures embedded in debris disks (Moro- Martin & Malhotra 2002), and the determination of 
characteristic wavelength ranges in which planetary radiation - in particular astrobiological 
tracers - can be separated from the underlying dust continuum radiation (Beichmann 1996, 
Fischer & Pfau 1997) are tasks which are beyond the recent study and will be considered in 
a subsequent investigation. 
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Fig. 1. — Absorption and scattering coefficients over the transition from Fe-deficient to Fe- 
rich Silicates. Arrow indicates increasing iron content from Mg Si0 3 — > Mg .95Feo.o5Si0 3 — > 
Mgo. 8 Feo.2Si0 3 -> Mg .7Feo.3Si0 3 -> Mgo. 6 Feo.4Si0 3 -> Mgo. 5 Feo. 5 Si0 3 -> Mg . 4 Feo. 6 Si0 3 -> 
MgFeSi0 4 — * Mgo.sFei.2 Si0 4 . The optical parameters have been calculated for grains with 
a radius of 1 /im. 
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Fig. 2. — Contributions from single dust compositions. Top: Mg Si0 3 (solid line), 
Mgo.95Fe .o5Si0 3 (dotted line), Mg .8Fe .2SiO 3 (dashed line), Mg a7 Fe .3SiO 3 (dash-dotted 
line), Mg . 6 Feo.4Si0 3 (dash-dot-dot-dotted line). Middle: Mg . 5 Feo.5Si0 3 (solid line), 
Mg .4Fe . 6 SiO 3 (dotted line), MgFeSi0 4 (dashed line), Mgo. 8 Fei. 2 Si04 (dash-dotted line), 
"Astronomical Silicate" (fat solid line). Crystalline Olivine (long-dashed line). Bottom: C 
(400 K - solid line), C (600 K - dotted line), C (800 K - dashed line), C (1000 K - dash-dotted 
line), Graphite (fat solid line). In all cases a single grain radius of a = 1/im and a disk mass 
of lO -10 ]^ are assumed. 
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Fig. 3. — Absorption and scattering coefficients for Mg-rich / Fe-deficient Silicate: MgSi0 3 . 

(a) upper row: grain radii 0.1/im (solid line), 1/im (dotted line), 10/im (dashed line); 

(b) lower row: grain radii 0.1mm (solid line), 1mm (dotted line), 10mm (dashed line). 
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Fig. 4. — Absorption and scattering coefficients for Mg-poor / Fe-rich Silicate: 
Mg . 8 Fe 1 . 2 SiO 4 . 

(a) upper row: grain radii 0.1/xm (solid line), 1/im (dotted line), 10/im (dashed line); 

(b) lower row: grain radii 0.1mm (solid line), 1mm (dotted line), 10mm (dashed line). 
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Fig. 5. — Absorption and scattering coefficients for Olivine (crystalline, but not oriented): 
Mgi. 9 Feo.iSi0 4 . 

(a) upper row: grain radii 0.1/xm (solid line), 1/im (dotted line), 10/im (dashed line); 

(b) lower row: grain radii 0.1mm (solid line), 1mm (dotted line), 10mm (dashed line). 
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Fig. 6. — Absorption and scattering coefficients for Carbon: 400 K configuration ("graphite 
poor"). 

(a) upper row: grain radii 0.1/xm (solid line), 1/im (dotted line), 10/im (dashed line); 

(b) lower row: grain radii 0.1mm (solid line), 1mm (dotted line), 10mm (dashed line). 
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Fig. 7. — Absorption and scattering coefficients for Carbon: 1000 K configuration ("graphite 
rich"). 

(a) upper row: grain radii 0.1/im (solid line), 1/zm (dotted line), 10/zm (dashed line); 

(b) lower row: grain radii 0.1mm (solid line), 1mm (dotted line), 10mm (dashed line). 
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Fig. 8. — Illustration of the relative contributions of light scattering and thermal equilibrium 
emission in the emergent SED. 

Solid line: Scattered + Reemitted radiation. Dashed line: Reemitted radiation only. 
Left column: Fe-rich Silicate (amorphous Olivine, MgFeSiO/t). 
Right column: Fe-poor Silicate (MgSiOs). 

Top: a = 0.1/im, Middle: a = 10/xm, Bottom: a = 1mm. Disk mass: 1O _1O M . 
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Fig. 9. — Dependence of the emergent SED on the disk density distribution. Disk density 
profile: q—1.0 (solid line), 1.5 (dotted), 2.0 (dashed), and 2.5 (dash-dotted). Top: a = 0.1/mi, 
Middle: a = 10/im, Bottom: a = 1mm. All calculations use amorphous olivine (MgFeSiO/t) 
composition and a disk mass of 1CT 12 M . 
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Fig. 10. — Influence of inner gaps on the SED. Inner disk radius: dust sublimation radius 

(solid line), 0.1 AU (dotted), 1AU (dashed), and 10 AU (dash-dotted). 

Left column: Fe-poor Silicate (MgSi0 3 ). 

Right column: Fe-rich (amorphous Olivine, MgFeSiO/t). 

Top: a = 0.1/zm, Middle: a = 10/im, Bottom: a = 1mm. Assuming a dust sublimation 
temperature of 1550 K, the sublimation radius of MgSi0 3 / MgFeSi0 4 amounts to 7.3 x 
10~ 3 AU / 1.3 x lO^AU (a = 0.1/im), 1.2 x HT 2 AU / 3.5 x 10~ 2 AU (a = 10/im), and 
3.3 x 10~ 2 AU / 3.2 x KT 2 AU (a = 1mm). Disk mass: 1Q- 1O M . 
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Fig. 11. — Influence of the outer disk radius on the grain size distribution. Outer radius of 
the disk: 25 AU (solid), 50 AU (dot), 100 AU (dash), 200 AU (dash-dot), 400 AU (dash-dot- 
dot-dot). Grain radius: 0.1/im (top), 10/zm (middle), 1mm (bottom). Amorphous olivine 
(MgFeSi0 4 ) composition and disk mass 1O _11 M are assumed. 
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Fig. 12. — Influence of the grain size distribution on the SED. 
Grain size distribution exponent: p = 2.5 (upper row) and 3.5 (lower row). 
Minimum grain size: a m i n = 0.1/im (left column), and 10/im (right column). 
Maximum grain size: a max — 0.1mm (solid line), 1mm (dotted), and 10mm (dashed). 
Amorphous olivine (MgFeSiO"4) composition and disk mass 10~ 10 M Q are assumed. 
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Fig. 13— The complete paper can be downloaded from 

http : //mc . caltech . edu/~swolf /downloads/ deb-sed . ps . gz. 

Distribution of normalized debris disk SEDs in the range [0.2/xm, 200/xm]. Ranges of vari- 
able disk parameters: p = [2.5,3.5], q = [1.0,2.5]. The particle size ranges from which 
the minimum and maximum grain size are randomly selected are given in each figure. 
Left /Right column: Minimum and maximum grain sizes are selected from a logarithmically 
equidistant / linear distribution of grain sizes. Color coding: Dark regions represent a high 
SED density, while light grey regions symbolize low SED densities. The black upper and 
lower curves in each figure mark the range within which solutions (SEDs) were found. 
See § 4.4 for a more detailed explanation and discussion. 
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Fig. 14. — SIRTF color-color diagrams corresponding to the SEDs shown in Fig. 2. [1] 
MgSi0 3 Mgo.gsFeo.osSiOs Mgo. 8 Fe . 2 Si0 3 Mg . 7 Fe .3SiO 3 - Mg . 6 Fe .4SiO 3 
Mgo.sFeo.sSiOg -> Mgo. 4 Peo.6Si0 3 ; [2] MgFeSi0 4 -> Mg . 8 Fe L 2SiO 4 ; [3] "Astronomical Sil- 
icate"; [4] Crystalline Olivine; [5] Carbon: [5.1] C (400 K), [5.2] C (600 K), [5.3] C (800 K), 
[5.4] C (1000 K); [6] Graphite. 
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Fig. 15. — SIRTF color-color diagrams corresponding to the SEDs shown in Fig. 9 (top 
left), Fig. 11 (top right), and Fig. 10 (bottom). Different line styles represent different dust 
grain radii: 0.1/zm (dotted), 10/zm (dot-dashed), 1mm (dashed). Top left: Different symbols 
represent different power-law exponents q of the disk density distribution: q — 1.0 (plus), 
q = 1.5 (asterisk), q = 2.0 (triangle), and q = 2.5 (square). Top right: Different symbols 
represent different outer disk radii: 25 AU (plus), 50 AU (asterisk), 100 AU (triangle), 200 AU 
(square), and 400 AU (rhombus). Bottom left/right: Different symbols represent different 
inner gap radii for the case of MgSi0 3 / MgFeSi0 4 grains: sublimation radius (plus), 0.1 AU 
(asterisk), 1.0 AU (triangle), and 10 AU (square). 
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Fig. 16. — SIRTF color-color diagrams corresponding to the SEDs shown in Fig. 12. 
Left/Right: Minimum grain radius amounts to 0.1/xm / 10/xm. Different line styles represent 
different maximum dust grain radii: 0.1mm (dotted), 1mm (dot-dashed), 10mm (dashed). 
Different symbols represent different power-law exponents p of the grain size distribution: 
p = 2.5 (plus), and q = 3.5 (asterisk). 



